function figout = zambia_figNumTr(sd_F1vec,mP,mPComp,mPN,mW,mPProf0,mProf0_C,sdf1)

%Graphs with baseline values of welfare parameters
set(0,'DefaultAxesFontName', 'Garamond')
set(0,'DefaultAxesFontSize', 16)

% Change default text fonts.
set(0,'DefaultTextFontname', 'Garamond')
set(0,'DefaultTextFontSize', 16)

%%2 4 1 3 %Do not need to divide column 2/4 by part because all part

figout = figure('NumberTitle','off','Position',[100, 100, 1200, 725]);
set(figout, 'PaperOrientation', 'landscape','Color', 'w')
subplot(2,2,1)
plot(sd_F1vec',mPN(:,1),'--ok',sd_F1vec',mPN(:,2),'k-'),...
    title('A. No Reward (R=0)'),axis([0 500 5 50]),ylabel({'Average follow-through','(trees per farmer)'})
    legend(strcat('Average follow-through: A=0 '),'Average follow-through: A=12 ','Location','northeast');
hold on
line([sdf1,sdf1],[5.01,50],'LineStyle',':','Color',[0.8,0.8,0.8])
hold off
sub_pos = get(gca,'position'); % get subplot axis position
set(gca,'position',sub_pos.*[1 1 1.09 1.05]) % stretch its width and height

subplot(2,2,2)
plot(sd_F1vec',mPN(:,3),'--ok',sd_F1vec',mPN(:,4),'k-'),...
    title('B. Large Reward (R=150)'),axis([0 500 5 50])
    legend('Average follow-through: A=0 ','Average follow-through: A=12 ','Location','northeast');
hold on
line([sdf1,sdf1],[5.01,50],'LineStyle',':','Color',[0.8,0.8,0.8])
hold off
sub_pos = get(gca,'position'); % get subplot axis position
set(gca,'position',sub_pos.*[1 1 1.09 1.05]) % stretch its width and height

subplot(2,2,3)
plot(sd_F1vec',mPN(:,1)./mP(:,1),'--ok',sd_F1vec',mPN(:,2)./mP(:,2),'k-'),...
    title('C. No Reward (R=0)'),axis([0 500 5 50]),xlabel('Uncertainty (standard deviation of F1)'),ylabel({'Average follow-through','conditional on take-up','(trees per farmer)'})
    legend(strcat('Follow-through given take-up: A=0 '),'Follow-through given take-up: A=12 ','Location','northeast');
hold on
line([sdf1,sdf1],[5.01,50],'LineStyle',':','Color',[0.8,0.8,0.8])
hold off
sub_pos = get(gca,'position'); % get subplot axis position
set(gca,'position',sub_pos.*[1 1 1.09 1.05]) % stretch its width and height

subplot(2,2,4)
plot(sd_F1vec',mPN(:,3)./mP(:,3),'--ok',sd_F1vec',mPN(:,4)./mP(:,4),'k-'),...
    title('D. Large Reward (R=150)'),axis([0 500 5 50]),xlabel('Uncertainty (standard deviation of F1)')
    legend('Follow-through given take-up: A=0 ','Follow-through given take-up: A=12 ','Location','northeast');
hold on
line([sdf1,sdf1],[5.01,50],'LineStyle',':','Color',[0.8,0.8,0.8])
hold off
sub_pos = get(gca,'position'); % get subplot axis position
set(gca,'position',sub_pos.*[1 1 1.09 1.05]) % stretch its width and height

end
